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Propagators in Coulomb gauge from SU(2) lattice gauge theory 



O 

o 

(N 



>: 

O 

o 

o. 

-(— > . 

Oh! 
(D . 



X 



Kurt Langfcld and Laurent Moyaerts 
Insitut fiir Theoretische Physik, Universitdt Tubingen 
D-72076 Tubingen, Germany. 
(Dated: June 15, 2004) 

A thorough study of 4-diniensional SU(2) Yang-Mills theory in Coulomb gauge is performed using 
large scale lattice simulations. The (equal-time) transverse gluon propagator, the ghost form factor 
d{p) and the Coulomb potential Vcoui{p) oc d?{p)f(p)/p^ are calculated. For large momenta p, the 
gluon propagator decreases like l/p^+'' with 77 = 0.5(1). At low momentum, the propagator is weakly 
momentum dependent. The small momentum behavior of the Coulomb potential is consistent with 
linear confinement. We find that the inequality Ocoui > o" comes close to be saturated. Finally, we 
provide evidence that the ghost form factor d(p) and f(p) acquire IR singularities, i.e., d{jp) oc 1/y'p 
and /(p) oc 1/p, respectively. It turns out that the combination g^doip) of the bare gauge coupling 
go and the bare ghost form factor do(p) is finite and therefore renormalization group invariant. 

PACS numbers: 11.15.Ha, 12. 38. Aw, 12.38.Gc 



I. INTRODUCTION 

Quark confinement is the key property which dictates 
the structure of matter at low temperatures and densi- 
ties. It is attributed to the low energy regime of Quantum 
Chromo Dynamics (QCD) and is, therefore, only accessi- 
ble by non-perturbative techniques. Lattice gauge theory 
provides a gauge invariant regularization of Yang-Mills 
theory, and corresponding numerical simulations do not 
require gauge fixing. 

More than twenty years ago, Mandelstam and 't Hooft 
pointed out that gauge fixing might be convenient for 
identifying the degrees of freedom which are relevant for 
confinement. Over the recent past, evidence has been 
accumulated by means of lattice simulations that topo- 
logical obstructions of the gauge field, such as monopoles 
and vortices, are responsible for confinement (see |lj for 
a recent review). Non-perturbative approaches based 
upon Dyson- Schwingcr equations (DSE) |^ IJ, varia- 
tional techniques ^, ^ ^ and flow equations 0, lal nec- 
essarily involve gauge fixing. These approaches address 
QCD Green functions which encode information on con- 
finement at low momenta. 

In Landau gauge, it was firstly put forward by Gri- 
bov Q and further elaborated by Zwanziger 10] that 
quark confinement is associated with a divergence of the 
ghost form factor in the infrared limit. First indications 
that the so-called horizon criterion for confinement is 
satisfied were reported in [ij. Using truncated DSEs, 
strong evidence for an IR divergent ghost form factor was 
reported in [Tj, [13 . Subsequently, many efforts were de- 
voted to detect the low energy behavior of the ghost form 
factor b y m eans of analytic JJj J^ J,.6, 17, 18, 21J and lat- 
tice [23,123,123 techniques. A good qualitative agreement 
between the DSE and lattice results was found. More- 
over, a tight relation between the vortex picture of con- 
finement and the Gribov-Zwanziger criterion was given 
in m. 



Coulomb gauge is most convenient for a variational ap- 
proach to Yang-Mills Green functions |j, |^, |^ . The cor- 
responding Hamiltonian contains an instantaneous inter- 
action between color sources, the so-called Coulomb po- 
tential Vcoui{r)- Hence, Coulomb gauge QCD offers the 
possibility to address the confining potential by pu rely 
analytic considerations. It has been shown recently |2g 
that the Coulomb potential is an upper bound of the 
static quark-antiquark potential. In pure Yang-Mills the- 
ory, both potentials are expected to be linearly rising at 
large distances r. At present, lattice simulations provide 
two different values for the asymptotic Coulomb force: A 
calculation of the potential involving the zeroth compo- 
nent of the gauge field yields acoui « (2 — 3) a [23, 123 , 
while acoui ~ cr |23| is reported if the Green function 
which defines the Coulomb potential is directly evalu- 
ated. 

The variational approach invokes a quasiparticle pic- 
ture for the (transverse) gluons. The trial wave func- 
tional generically is Gaussian, and the gluon dynamics is 
encoded in the (quasi particle) gap function luqp (which 
is derived in a gluonic quasiparticle picture from the in- 
verse of the instantaneous gluon propagator). It was re- 
cently proposed to supplement the Faddeev-Popov de- 
terminant (see below) to the wave functional in order 
to strengthen the impact of configurations close to the 
Gribov horizon igj] . The latter modification of the wave 
functional has certainly a strong impact on the gap func- 
tion at least at small momenta: while it is suggested in J^ 
that ujQp approaches a constant in the IR limit, an IR 
divergence was reported in [g. Both approaches employ 
a truncation of the resulting Dyson equations which does 
not account for wave functional renormalization. 

In the present paper, we perform a thorough lattice in- 
vestigation of the transverse gluon propagator, the ghost 
form factor and the Coulomb potential. Our results sug- 
gest that (at least in four dimensions) the equal time 
gluon propagator is weakly momentum dependent in the 
IR regime. For the first time, wave function renor- 



malization constants for the gluon and ghost fields are 
obtained. Our lattice results indicate quite a sizable 
anomalous dimension for the gluon propagator yielding 
Gdpl) — > I/IpT"'''' where 77 « 0.5. Our results for the 
Coulomb potential are compatible with linear confine- 
ment and suggest that the inequality acoui > cr is almost 
saturated. This finding favors the result in [23. We will 
find, however, that the result acoui ~ (2 — 3) a, reported 
in |23,Ell, cannot be ruled out. 



where A° are the gauge fields, F° denotes the usual field 
strength tensor and J^ is an external source. 
Det[— V ■ D] is the Faddeev-Popov determinant (D is the 
gauge covariant derivative), and go the bare gauge cou- 
pling. Equation |^ will serve as the starting point for 
the lattice simulations reported in this work. 

The canonical momenta conjugated to the gauge fields 
are introduced by 



II. 



THE PARTITION FUNCTION 



Let us briefly review the relation between the func- 
tional integral formulation and the Hamilton formulation 
of Coulomb gauge Yang-Mills theory following ^J . Our 
starting point will be the generating functional for Eu- 
clidean Green functions in Coulomb gauge, i.e., 

Z[J] = /"pAexpj /"/:dx|5(5,AJ.)Det[-V-D],(l) 



VU exp 



dx ( -n^ 
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exp -- / dxFg^F'^ 



Eliminating the non-dynamical part of the gauge fields, 
Aq, implements Gauss's law in the partition function. 



Z = VA VU exp 



IlfA' 



i(tf 



B^) - s„A° ■ J"] I s(d,At) i([D.,nT - SoJS) Det[-V - D]. (2) 



Let us introduce the transverse and longitudinal parts of 
the canonical momenta by 11'' — U"^ — S/cj)"'. The inte- 
gration measure factorizes, i.e., VU ~ VU^Vcl), and the 
integration over (j) cancels the Faddeev-Popov determi- 
nant. One finally obtains 



(3) 



Z = I?A_LPni exp <^ / {i n'i^^A'-'' - H 
Thereby, the Hamiltonian Ti± is given by 



Hj 



d-'x 



i(ni 



B^ 



5oAj 



1 
2C^ 



d^x d^y p(x) Vcoui(x, y) p(y), (4) 



where C2 is the quadratic Casimir of the gauge group 
SU(2), i.e., C2 = 3/4, and B is the chromo-magnetic 
field. We have introduced the color charge density by 

r''^A'i{x)U'i^{x) + Jo{x). 



p{x) 



(5) 



The "tree level" term which mediates the interaction be- 
tween color charges, i.e., 

Vcoui(x,y) = -C2 gl A'f-i(-A)M-i|(,,y) , (6) 

will give rise to the so-called Coulomb potential upon 
averaging over the gauge fields. Thereby, 

M = -V-D (7) 



is the Faddeev-Popov matrix. 



III. THE NUMERICAL APPROACH 

A. Setup 

Configurations of a N'^ cubic lattice with lattice spac- 
ing a are generated using the standard Wilson action. 
A lattice configuration is represented by a set of unitary 
matrices U^{x) e SU{2). The size of the lattice spacing 
as a function of the Wilson /3 parameter is obtained from 
the interpolating formula p4| 
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where 



/3o - y 



/9i = ^ 



(9) 
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The first two terms on the r.h.s. of © are in accordance 
with 2-loop perturbation theory. The term d/ (3 repre- 
sents higher-order effects and the term c is a dimension- 
less scale factor to the string tension. It was observed |2J| 
that the choice 



4.38(9) , 



1.66(4) 



(10) 



reproduces the measured value acr to very good preci- 
sion. 

Two possible definitions of the gauge field ^t(a;) in 
terms of the link matrices are explored in the present 
paper. Decomposing a particular SU(2) matrix by 



U^{x) = ul{x) + ^<(x)r^ 



(11) 



where r" are the Pauli matrices, the standard definition 
of the gauge field is given by 



^go4(x) + 0{c?) = 2ul{x) 



(12) 



Noticing that the continuum gauge fields actually trans- 
form under the adjoint representation, the definition 

ag^Al{x) + 0{c?) = 1ul{x)ul(x) (13) 

accounts for this particular transformation property |3lj . 
Note that both definitions coincide in the continuum 
limit a — > 0. In the context of minimal Landau gauge, 
both definitions also coincided at the practical level 
within the scaling window |24|. 

Finally, we define the momentum on the lattice for the 
particular direction /i 



2 /TT 



(14) 



where —N/2 < n^ < N/2 labels the Matsubara fre- 
quency. This definition minimizes rotational symmetry 
breaking for momenta close to the boundary of the Bril- 
louin zone. 



B. Propagators and Coulomb potential 

Once and for all, we choose a given time slice of the 
four dimensional space-time by fixing t = to, and consider 
propagators which are defined within the emerging 3- 
dimensional hypercube. 

The (bare) gluon propagator is defined by 

G°/^x-y) := (A^(x,io)4(y,t„)). (15) 

The corresponding propagator in momentum space, i.e., 

G°/'(p)-/rf'a,-G°/''(x)e*P-, 

is transversal by virtue of the gauge condition 
9iyl^(x, to) = and diagonal in color space, i.e., 



G^-(p) = (%-^) 5^'G\p), (16) 



'V p2 



G'O(p) = 4^ 



(17) 



The dimensionless quantity /° is the (bare) gluon form 
factor. Since the gauge potential has energy dimension 
one, the quantity G°(p) has the dimension of an inverse 
energy. 

The (bare) ghost propagator is defined as the expecta- 
tion value of the inverse Faddeev-Popov operator M in 
0, i.e.. 



D°"^(x-y) = (M-1[A]P'' 



(x,y) 



Since D is diagonal in color space, we may write 



(18) 



D°''\p)^6''''D°{p), D°{p)^ 



d°(p) 



in momentum space. Thereby, we have introduced the 
(bare) ghost form factor dP{p). Since the quantity 
(A/^^[A]|(x,y)) has energy dimension 1, the ghost form 
factor dP{p) is dimensionless and reduces to unity for a 
free theory. 

The Coulomb potential is given by the expectation 
value of the expression ©: 



-C, gl/[AI-'[A]i-A)M-'[A]]1l^ . (19) 



1 ab 



For the Fourier transform of the Coulomb potential, we 
make the ansatz 



^^CoullPJ - ~G2 5o f— T^ 



(20) 



where dP{p) is the bare ghost form factor. The dimen- 
sionless function /(p) measures the deviation from the 
factorization 



M-^[A] (-A)A/-i[yl] 

{m-'[A]) i~A) {m-'[A]) , (21) 
in which case /(p) ~ 1. 



C. Renormalization 

Renormalization of Yang- Mills theories in four dimen- 
sions implies that the bare coupling acquires a depen- 
dence on the ultraviolet (UV) cutoff Auv, i.e.. 



go -^ goi^UV / ^scale) ■ 



(22) 



Rather than the bare coupling constant, the Yang-Mills 
scale parameter Ascaie plays the role of the only param- 
eter of the theory. In the context of (quenched) lattice 
gauge simulations the string tension a is widely used as 
the generic low-energy scale. In this case, the cutoff de- 
pendence of the bare coupling is implicitly given by the 



/3 dependence of a a? {[3) where j3 = 4/(7q is related to 
the cutoff by A^/v — 7T/a{(3). Finally, the relation be- 
tween aa^ and /3 is provided by the formula ^ which 
interpolates the calculated values. 

In addition, wave-function renormalization constants 
generically develop a dependence on Auvf-^scaie- The 
lattice bare form factors of the previous subsections, fg 
and dP , depend on the momentum p'^ and on the UV cut- 
off Auv (given in units of the string tension) or, equiva- 
lently, on the lattice coupling f3. The renormalized form 
factors are obtained upon multiplicative renormalization 



rrip',n') = z-\p,^i)fg{p\p) (23) 

d-"(/,/i2) = Z^\p,^Ji)d\p\p) (24) 

using the renormalization conditions 

/-"(M^A^^) = 1, d-V,M') = l- (25) 

Finally note that the Coulomb potential Vcoui(p) H2U|) 
appears as a convolution of propagators. Let us factor 
out the wave function renormalization of the ghost fields, 
and let us introduce the renormalization constant Zf of 
the 4-point vertex function involving four ghost fields: 

^Coui(PW^') = ^7\/3,m)^3"'(Am)^^c°oui(p,/3), (26) 

where V^^^^ i^ given in (|19|) . If the theory is renormal- 
izable without a 4-ghost counter term at tree level, then 



lim Zf{f3,fi) = constant 



(27) 



and the "factorization" function /(p), which is implicitly 
defined by l(^ . does not acquire an UV divergence. Be- 
low, we will suppress the superscript "ren" and only deal 
with renormalized quantities. 

Multiplicative renormalizability of the theory implies 
that a rescaling of the data for each /3 value (indepen- 
dently of the momentum) is sufficient to let the form 
factors fall on top of a single curve describing the mo- 
mentum dependence of the corresponding renormalized 
quantity. In practice, suitable "matching factors" are de- 
termined which "collapse" data obtained at different (3 
on a single curve. The matching factors are then directly 
related to the wave function renormalization constants. 
This procedure is described in detail in (3^ Sec. V.B.2] 
and HJI. 



D. Gauge fixing 

Let us consider the time slice ip in which we define the 
transverse (equal time) gluon propagator H15|l and the 
ghost propagator, respectively. 



The lattice configurations are generated without any 
preference for a particular gauge. Subsequently, the 
gauge transformations n{x) are adjusted according to 



Fu[n] = 



where 



x,i—1...3 



RcTr l-[/^^(f,to) 



min, (28) 



u^{x) ^ n{x) Uf,{x) n^x + n) . 

It is well known that the condition (|28|l fixes the gauge 
up to gauge transformations which only depend on time, 
i.e., n{t). Note, however, that the average over the un- 
fixed gauge degree of freedom does not affect the (equal 
time) gluon propagator. The same is true for the ghost 
propagator. 

The calculation of the non-instantaneous gluon propa- 
gator, i.e.. 



G: 



gen ab 



(x-y,i-io) := (A^(x,i)A5(y,to)) (29) 



would inevitably require complete gauge fixing, since the 
residual gauge freedom would imply 



Gf™'^''(x-y,i-io) = 



for t ^ to 



The standard procedure is to involve the gauge field 
Aq{x) for residual gauge fixing. An example is the 
Cucchieri-Zwanziger condition [33J 



Ft[n] = ^ReTr(l-C/^^(f,i) 



n{t) 



(30) 



x,t 



It is also possible to fix the residual gauge freedom with- 
out incorporating Aq{x). In this case, the Aq{x) inte- 
gration in section ^] can still be performed thereby en- 
forcing Gauss's law (see Q)- The advantage is that the 
tight relation between the functional approach (Q and 
the Hamiltonian formulation (Q is preserved. 

After averaging over the gauge fields, the Fourier trans- 
form of 129|) can be written as 

Gff-\p,po) = 5^' (<5,, - ^) ^^ , (31) 

where the dimensionless quantity i^(p,Po) is a form fac- 
tor. If the residual gauge degree of freedom is not fixed, 
the gauge fields corresponding to different time slices 
would be uncorrelated, and the form factor would not 
depend on pQ. We expect that the residual gauge fixing 
only introduces weak correlations between gauge fields 
of different time slices implying that the form factor only 
slightly depends on po- A detailed lattice investigation 
is left to the near future. In the present paper, we will 
focus on the equal time gluon propagator (|15|l . 

In practice, we used a simulated annealing method to 
locate the minimum of H28|l . The idea is to consider the 



gauge functional F[7[ri] as the action of a field theory 
with respect to the set of the gauge transformations ri(x), 
whose partition function is given by 

ZsA = J nneicpi-PsAFum). 

By increasing step by step the free parameter /3s A (which 
plays the role of the inverse temperature), one tries by 
cooling to retrieve the ground state of the fictitious field 
theory. This procedure corresponds to the minimization 
of Fu [g] . The thermalization steps can be performed us- 
ing the standard Creutz update algorithm supplemented 
by microcanonical reflection steps. After the domain of 
attraction has been located by simulated annealing, we 
perform iterated overrelaxation to ensure the desired pre- 
cision of the gauge condition. Details will be presented 
elsewhere. At least in minimal Landau gauge, it has 
turned out that the present method of gauge fixing is 
robust against the impact of Gribov ambiguities [2J| . 



IV. NUMERICAL RESULTS 



propagator, they concluded that their data, obtained for 
momenta below 2 GeV, did not reach the perturbative 
regime. In the present paper, the high momentum regime 
is also explored. It turns out that the trend with 77 as 
large as 0.5 continues at least up to momenta as large as 
12 GeV. 

Also shown in figure ^ is the cutoff dependence of the 
wave function rcnormalization constant Z3. Unfortu- 
nately, the achieved numerical precision does not allow 
for a definite conclusion on the functional form which 
parameterizes the cutoff dependence of Z3 . 

In the context of a variational approach Q, IM 0' 
1/G(p) is interpreted as energy dispersion relation of con- 
stituent gluon fields. Note, however, that the full gluon 
propagator G(p) which comprises all non-pcrturbative 
effects is parameterized by G(p) — g(p)/'^Qp(p), where 
(7(p) accounts for effects which are beyond the leading or- 
der of the quasi-particle (QP) picture. Hence, the inter- 
pretation of the inverse propagator as dispersion relation 
must be taken with care. Our results for the inverse gluon 
propagator 1/G(p) for several lattice sizes are shown in 
figure [3 



Gluon propagator 



Lattice simulations were carried out using 26'*, 
32^ and 42"* lattices. Inverse couplings /3 G 

{2.2, 2.3, 2.4, 2.5, 2.6, 2.7, 2.8} were employed. Within the 
momentum window, very good scaling is observed imply- 
ing that cutoff effects (due to the finite value a) and finite 
size effects are small. Our final result for the transverse 
equal-time gluon propagator G(p) H17|l is shown in figure 
n Only data for the 42'* lattice are shown for clarity. We 
have checked that the data for the 26"', 32"* lattices fall on 
top of the same curve. As a renormalization condition, 
we have chosen 



G(p = lGcV) = f [GeV]" 



(32) 



for a renormalization point ti — 1 GeV. For comparison, 
we have also shown the data obtained by Cucchieri and 
Zwanziger in [33 ■ These data are obtained for j3 = 2.2 
and 28^, 30* lattices by using an iterated overrelaxation 
method for gauge fixing. We find good agreement. 

At small momentum, the propagator becomes roughly 
momentum independent and seems to approach a con- 
stant in the IR limit |p| — > 0. At large momentum, the 
momentum dependence is well approximated by 



G(p)« 



1 



^qcd \ 
IpI ) 



7? = 0.5(1) 



(33) 



where Aqcd is the renormalization group invariant scale 
parameter. The solid line in figure Q] corresponds to a fit 
with rj = 0.5. Also shown are the lattice data from [S^l 
(by courtesy of Attilio Cucchieri). The large value 77 
was already anticipated in [33 ■ Since the authors fo- 
cused their investigations on the IR behavior of the gluon 



B. Ghost form factor 

In order to obtain the ghost propagator IjlHII . an in- 
version of the Faddeev-Popov matrix is performed after 
(possible) zero modes have been removed jlj ■ This pro- 
cedure requires the solution of large, but sparse, systems 
of equations. In order to remove zero mode contributions, 
we solve 

MMx = Mb 

rather than the system Mx = b (see the method proposed 
in dlij). For this purpose, we use the minres algorithm 
with Jacobi preconditioning |3J| . Our final result for the 
ghost form factor d(p) obtained from simulations using 
a 26* lattice and /? G {2.15,2.2,2.3,2.4,2.5} is shown in 
figure m right panel. We observe perfect scaling. To be 
specific, we choose 

d(p==3GeV) = 1 

as a renormalization condition. 

In order to explore the high momentum dependence, 
we make a logarithmic ansatz supplemented with an 
anomalous dimension 7g, i.e.. 



d(p) = 



ln(^|p|/AQCDJ 



P > Aqcd • (34) 



The UV fit shown in figure [5] (right panel) corresponds 
to 



1.03(2) 



7go 



0.26(2) 
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FIG. 1: The transverse equal-time gluon propagator as a function of the momentum: our data (A'^ — 42) and data from |3^ 
(left panel). The gluon wave function renormalization constant Zs as a function of the UV cutoff. Also shown are the wave 
functional renormalization constant of <7o'^(p) and Zf 12611 : the y-axis is arbitrarily scaled (right panel). 
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FIG. 2: Momentum dependence of the inverse gluon propagator 1/G'(p) (left panel). Ghost form factor (right panel). 



Since the fitting parameters are strongly correlated, it 
is difficult to extract a reliable value for Aqcd- Taking 
also into account the high momentum behavior of the 
Coulomb potential (see below), we find that 



Aqcd = 0.96(5) GeV 



reproduces the UV data. Such a high value for the Yang- 
Mills scale parameter seems to be generic for the lattice 
regularization (see pH). For the IR analysis, we adopt a 



simple scaling law: 
rf(p) = 



p|VA2 



QCD 



P < Aqcd • (35) 



The IR fit is also shown in figure |5] (right panel). We find 



a,r ~ 1.55(1) , 



K = 0.245(5) 



(36) 



Hence, our lattice results suggest that the ghost form 
factor roughly diverges like l/-\/|p| in the IR limit |p| ^ 
0. 




2 3 

p [GeV] 



2.5 



oo 



1.5 



> 



0.5- 



- 




1 ' 


' 1 ' 1 


o 


b=2.5 

b=2.4 

b=2.3 

b=2.2 

b=2.15 

peit. theory 






- 


i 


*«?.««^' 


^ - 


^ 




1 


1 1 1 1 



2 3 

p [GeV] 



FIG. 3: p^ Vcoui(p) (Coulomb potential) as a function of the momentum p (left panel), p^ Vcoui(p)/87r in units of the string 
tension a (right panel). 



Let us explore the ghost wave function renormahza- 
tion constant Z^. If one apphes the momentum match- 
ing technique to ^g '^o(p, P), one finds that the matching 
factors are constant within the numerical precision, see 
figure n right panel. Our findings therefore suggest that 
the bare coupling squared times the bare ghost propaga- 
tor, i.e., 

gl D^o'iP) , (37) 

is a renormalization group invariant. 



C. Coulomb potential 

In the present approach, the Coulomb potential ap- 
pears as a convolution of Green functions (see H19|l ). 
In full Yang-Mills theory, wave functional (and vertex) 
renormalization applies (see H2f)|l 'l. and the prefactor of 
the potential must be specified by a renormalization con- 
dition. The prefactor can be determined by demanding 
that the perturbative result is recovered [S^ at large mo- 
mentum, i.e., 



P^ Vcoui(p) 



677 



lllnlplVA^P^ 



(38) 



for IpI ^ ^QCD- An elegant method which circumvents 
the cumbersome determination of the prefactor was put 
forward in EHH- 



Our numerical findings for p^ Vcoui(p) are summarized 
in figure 131 (left panel). At large momentum, the lattice 
data nicely show the logarithmic correction. The data 



of low energy regime are well reproduced by the scaling 
ansatz 

c 



p Vcoul(p) 



(39) 



where c is related to the Coulomb string tension acoui 
by (Tcoui = cAq^^/Stt in case Vcoui is compatible with 
linear confinement {6 = 1). Here, we find (see figure |21 
left panel) 

6 = 1.05(5) , c = 4.9(2) . (40) 

Also shown in figure [21 (right panel) is the combination 

p" Vcoui(p) /Sttct 

which should approach ctcoui/c in the limit |p| -^ 0. Al- 
ready for p > 2 GeV, the perturbative result is recovered 
to good precision. A plateau is reached at the inter- 
mediate momentum range 1 GeV < p < 2 GeV. The 
function is slightly increasing again for p < 1 GeV. The 
data within the observed momentum window are sta- 
ble against finite size effects. An extrapolation p — > 
is cumbersome. Our findings suggest that the inequal- 
ity a < CTcoui is almost saturated. This conclusion is in 
agreement with the results published in [23 ■ However, 
we point out that values Ucoui/o' ranging from 2 to 3 re- 
ported in [23, I23I cannot be ruled out from the present 
data. Larger lattices and higher statistics will be neces- 
sary to explore the deep infrared regime. 

D. Factorization 

From the results of the previous subsections it is al- 
ready clear that the function /(p) cannot weakly depend 



with 




./ 



2 3 

p [GeV] 



2.7(1) 



Kf = 0.58(5) 



(43) 



Our findings suggest the following interpretation: due 
to asymptotic freedom, the expectation value on the 
left hand side of H21|l factorizes in the high momentum 
regime. At low momentum, non-Gaussian correlations 
between the ghost fields play an important role and the 
factorization assumption is ruled out. 

Finally, we point out that the sum rule, dictated by 
perturbation theory, i.e., 2jgo + 7/ = 1, is satisfied to 
good precision: 



2 7go+7/ = 0.99(12) 



(44) 



This serves as a consistency check of the lattice renor- 
malization procedure. 



FIG. 4; The factorization function /(p). 



on |p|. Factorization such as H21|l is generic if the fields 
arc distributed according to a Gaussian probability dis- 
tribution. Since the ghost fields are strongly interacting 
at low momenta, one expects that /(p) significantly de- 
pends on the momentum in the IR regime. 

In order to extract /(p) from the lattice data we ap- 
plied the following procedure: Dividing the bare data for 
P^^coui ^ith the bare ghost form factor squared provides 
information on the unrenormalized function .f°(p). Ap- 
plying the momentum matching technique 2A\ yields the 
renormalization constant Zf and the renormalized func- 
tion /(p). Our numerical result for Zf as a function of 
the UV cutoff is consistent with (see figure^ right panel) 

lim Zf(/3,fi) = constant. 

The function /(p) is shown in figure 01 We observe a 
weak momentum dependence in the UV regime, and an 
IR singularity at small momentum. The high energy data 
are reproduced by the ansatz 



/(P) 



7/ 



P > Aqcd ■ (41) 



ln(|p|/AQCD 

The UV fit shown in figure 0] corresponds to 

aL « 1.05(4) , 7/ = 0.47(3) . 

The low momentum data support the existence of a l/|p 
singularity, i.e.. 



/(P) 



./ 



(IpP/a 



QCD 



p < Aqcd ■ (42) 



CRITICAL REMARKS 



Equal time gluon propagator 



In subsection IIV Al the lattice data for the transverse 
equal time gluon propagator G(p) are consistent with the 
high momentum behavior G'(p) oc |p|~^^'', r] — 0.5(1). It 
is tempting to conclude that asymptotic freedom implies 
the decrease of the equal-time gluon propagator G(p) like 
l/|p| for large momenta. We stress that this conclusion 
does not necessarily apply: The equal-time gluon prop- 
agator G(p) is derived from the generalized propagator 
(|31|l by momentum projection i.e.. 



G(p) ex Z-' JdpoGS°''ip,po) 
F{p,Po) 



= Zr, 



dpo 



(45) 



After proper renormalization, the dimensionless form fac- 
tor can be parameterized by 



F{p,Po) 



Aqcd Pa 



(46) 



where Aqcd is the low energy scale parameter (e.g., 
given by the string tension). For large momenta |p| S> 
Aqcd, one obtains 



G(P) 



1 



du /i(0, w) . 



(47) 



If the function h{x,y) is completely regular and bounded 
for x,y E [0, cxd[, the integral in the latter equation would 
exist implying that G(p) oc t^ for large |p|. However, 
the form factor _F'(p,po) weakly depends on po by the 
choice of the gauge (see discussions below (|31f) '). Let us 



speculate that for this reason we might assume for u ^ 1 
that 



h(0, u) oc — 7- 



1 



u^-'' 



< 77< 1 



In this case, the asymptotic form of the equal time prop- 
agator would be given by 



G(p) (X 



1 /A 



IpI J 



where the divergent factor {^uv / AqcdV has been ab- 
sorbed by the wave function renormalization constant Z^. 
The above picture is only one scenario which would ex- 
plain the lattice data. An lattice investigation of the 
po, P dependence of the generalized gluon propagator 
Gf"'° (p,Po) in (|31() will reveal whether the above sce- 
nario applies for Coulomb gauged YM-theory. Such an 
investigation is left to future work. 



B. Renormalization group 

Perturbation theory in Coulomb gauge j35l l3q is 
plagued by superficial divergences originating from in- 
stantaneous loops. When "supplementary rules" are 
introduced to define zero momentum integrals, a self- 
consistent framework is obtained. Within this approach, 
the Zg and the ghost renormalization constants, Zc and 



Z(j, are related by [3£ 



Zq Zc — 1 , 



ZlZc 



Zc = 



1 - 



7.9^ 



127r2 



(48) 



We point out that our lattice results indicate that the 
product of bare gauge coupling 170 and bare ghost form 
factor Z?07 i-e., gi Dq, is finite, which contradicts (|48|l . 
Further studies, involving simulation with larger lattices 
on one hand and partial resummations of perturbation 
theory to handle instantaneous loops on the other hand, 
seem necessary to resolve the above discrepancy. 



Let us focus here on the high and the low momentum 
behavior. For |p| ^ 1 GeV, our results for the (equal 
time) gluon propagator show a large anomalous dimen- 
sion, i.e., 



G(p)(xl/|p|H 



?7« 0.5(1) . 



The high momentum behavior of the ghost form factor 
and the factorization function /(p) show the character- 
istic logarithmic momentum dependence, i.e., 



d(p) oc In 



/(P) ex In 



A, 



QCD 



A, 



QCD 



-11 



-Igo = 0.26(2) 
If = 0.47(3) , 



while a scale parameter of Aqcd — 0.96(5) GeV is con- 
sistent with the lattice data. The high momentum behav- 
ior of the Coulomb potential known from perturbation 
theory |33, i.e.. 



P^ Vcoui(p) 



Gn 



llln|p|VA|^^ 



is recovered to high precision. The anomalous dimensions 
extracted from the lattice data are consistent with the 
sum rule 2 jgo + 7/ = 1- 

In the IR limit |p| -^ 0, the quantities d(p) and /(p) 
develop singularities. The data for the Coulomb potential 
are parameterized in the IR regime by 



Vcoul(p) 



Stt acoui 



and are consistent with linear confinement. We find that 
the inequality 



VI. CONCLUSIONS 

A thorough lattice study of the equal-time propagators 
of four dimensional SU(2) Yang-Mills theory in Coulomb 
gauge has been performed. Our simulations extend pre- 
vious studies to larger lattices and number of different 
lattice spacings. Our lattice data are in good agreement 
with those of earlier publications [23, 123, |33| ■ 



We studied the transverse (equal time) gluon propaga- 
tor G(p), the ghost form factor d{p) and the Coulomb 
potential Vcoui(p) H20|l . With the help of the momen- 
tum matching technique [2J|, the wave functional renor- 
malization constants are obtained for the gluon and 
ghost fields. In particular, we find that the combina- 
tion (7q (io(p) of the bare gauge coupling go and the bare 
ghost form factor do is renormalization group invariant. 



is almost saturated in agreement with [23 • We stress, 
however, that values in the range (TcouI = (2 — 3) cr as 
reported in |27l l2q are not ruled out by our data. The 
singularity of Vcoui(p) « l/|p|^, responsible for linear 
confinement, arises from IR singularities of the functions 
d(p) and /(p): 



d{p) ex 



1 



,10.49(1) 



/(P) « 



1 



11.17(10) 



The complete momentum dependence for the above func- 
tions can be found in section HVI 
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